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Understanding the nature of the interaction at the graphene/metal interfaces is the basis for graphene-based 
electron- and spin-transport devices. Here we investigate the hybridization between graphene- and 
metal-derived electronic states by studying the changes induced through intercalation of a pseudomorphic 
monolayer of Cu in between graphene and Ir(l 11), using scanning tunnelling microscopy and photoelectron 
spectroscopy in combination with density functional theory calculations. We observe the modifications in 
the band structure by the intercalation process and its concomitant changes in the charge distribution at the 
interface. Through a state-selective analysis of band hybridization, we are able to determine their 
contributions to the valence band of graphene giving rise to the gap opening. Our methodology reveals the 
mechanisms that are responsible for the modification of the electronic structure of graphene at the Dirac 
point, and permits to predict the electronic structure of other graphene-metal interfaces. 

The discovery of the unique transport properties of graphene 1,2 has stimulated the search for practical 
applications of this material in devices based on the transport of electrical charge or/and spin. In most of 
these devices, the graphene-metal interfaces are utilized as a contact for charge or spin injection. Hence the 
crystallographic and electronic properties of such interface determine the charge/spin injection efficiency of the 
corresponding contact, which is one of the main motivations for studies of the graphene-metal interface 3 5 . 

From another, more fundamental point of view, the study of the bonding mechanism at the graphene-metal 
interface is a very interesting problem in itself, and up to now such systems are far from the being fully 
understood 6 . Several factors influence the electronic properties of the graphene-metal interface: charge transfer 
from/onto graphene-derived n states, hybridization of the electronic valence band states of graphene and the 
metal, and the lattice match between the graphene and metal surface. These factors determine the behaviour of the 
n states in the vicinity of the Fermi level, E F , (as, e. g., a deviation from the linear behaviour characteristic for free- 
standing graphene) as well as the appearance of an energy gap in the spectrum of the graphene-derived electronic 
states at the so-called Dirac point, E D . From the point of view of electronic structure, two scenarios may be 
distinguished 6 : those in which the characteristic linear dispersion near E D is largely preserved, such as on Ir(l 1 1), 
Pt(lll), and Cu(lll), and those where a massive rearrangement of bands occurs, such as Ni(lll), Co(0001), 
Ru(0001), for example. 

Intercalation of metals in between graphene and substrates offers an interesting scientific playground to 
investigate the metal-graphene interaction. First, graphene may become decoupled from strongly interacting 
substrates, such as in the case of Au or Al intercalation in between graphene and nickel 710 . Such intercalated 
layers may also change the carrier concentration in graphene, and even change the carrier type (from electrons to 
holes) such as in the case of Au intercalation in graphene/SiC(0001) u . Intercalated metals may also enhance the 
magnetic coupling between a ferromagnetic substrate and graphene, with a view to utilizing graphene as a spin 
filter 12,13 . Moreover, the passivating, protecting function of graphene may also be used in such systems 14,15 . The 
intercalated layer in itself may bring new properties to graphene, such as in the case of lithium where supercon- 
ductivity has been predicted to occur 16 . 

The different factors that determine the bonding between graphene and metal can be analyzed by studying 
graphene-based intercalation systems, and recently, such attempts were undertaken for graphene on Pt(l 1 1) and 
Ir(lll), and on pseudomorphic layers of Co and Ni grown on these substrates 1719 . On top of Pt(lll) or Ir(lll), 
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graphene has properties almost like the free-standing phase, as 
judged by its electronic structure derived from photoemission and 
scanning tunnelling spectroscopy 20,21 . Intercalation of 1 ML of Co or 
Ni leads to a strong buckling of the graphene layer similar to the one 
formed on Ru(0001) 22 , and a large energy gap between n and n* 
states around the K point occurs, due to the broken symmetry for 
the two carbon sublattices in the graphene unit cell, induced by the 
strong hybridization of the graphene n and Co, Ni 3d valence band 
states; in both cases the linear dispersion of the graphene -derived 
states in the vicinity of E F is not conserved 1719 . 

A very different situation is found for intercalation of noble metals 
(Cu, Ag, Au) and the formation of single close-packed pseudo- 
morphic layers. Here, due to the absence of d states in the close 
vicinity of E P and the change of the doping of graphene upon inter- 
calation, the influence of n — d hybridization effects on the electronic 
structure of graphene is much weaker 7,8,23,24 , and the linear dispersion 
of the graphene-derived n states survives. The intercalation of single 
layers of noble metals is a suitable model system to investigate the 
interaction between graphene and metals, since it permits to follow 
the competition between the different (substrate and intercalated 
layer) electronic states. 

Here we study a pseudomorphic single layer of Cu atoms inter- 
calated in between graphene and Ir(lll), using angle-resolved 
photoelectron spectroscopy (ARPES), scanning tunnelling micro- 
scopy (STM), and state-of-the-art density functional theory (DFT) 
calculations. The interaction between the metal and graphene can be 
studied in detail because of the sharp and well separated signatures of 
the Cu 3d bands in photoemission. Our results yield a complete 
picture of hybridization between the Cu and carbon derived states 
and the consequent opening of a band gap in the graphene n bands; 
they are of importance for the understanding of the bonding mech- 
anism at graphene-metal interfaces in general. This system is also 
interesting since it permits the creation of a single layer of Cu under 
considerable tensile strain, because of the lattice mismatch of 6.2% 
between the two metals. If deposited on top of an Ir (1 1 1) surface, this 
strain leads to de-wetting and the creation of 3D islands/clusters even 
at low coverages, or to dendritic growth of layers as in the case of Au 
on Ir(lll) 25 or incomplete growth of Cu on Pt(lll) 26 . 

Results 

Graphene layers on Ir(lll) were prepared in a standard way via 
cracking of propylene gas at 1100°C STM images demonstrate the 
long range ordering of the moire structure over several hundreds nm 
as shown in Fig. 1(a); this image is an extract of the originally 
acquired 300 X 300 nm 2 STM data set. The lower inset shows an 
atomically resolved STM image of the moire unit cell of graphene/ 
Ir(l 1 1), which demonstrates the periodicity of 10 graphene unit cells 
over 9 unit cells of Ir(lll), consistent with the corresponding low 
energy electron diffraction (LEED) pattern shown in the upper inset 
of Fig. 1(a) 27 29 . At the bias voltages used in the experiment (Ut = 
+ 0.3 V and Uj — +0.5 V), graphene/Ir(lll) is imaged in the so- 
called inverted contrast 27,30 , where the crystallographically highest 
ATOP positions are imaged as dark areas and the lowest HCP and 
FCC places are imaged as bright ones, demonstrating the strong 
influence of the electronic structure on STM contrast (for the defini- 
tions of the high-symmetry positions of the graphene/metal moire 
structure see Fig. SI in the supplementary material). 

Intercalation of a 1 ML-thick Cu layer underneath graphene was 
performed via stepwise annealing of a deposited layer of Cu (with a 
nominal thickness of 1.5 ML) on graphene/Ir(lll). During this pro- 
cedure, the intensity of the "live" C Is and Ir 4/photoelectron spectra 
was taken as the intercalation proceeded as a function of annealing 
temperature [see Fig. S2(a,b) in the supplementary material]. This 
method permits to follow the formation of an intercalated Cu layer 
underneath graphene on Ir(lll), and the absence of any additional 



low binding energy (BE) components in the C Is spectra after inter- 
calation indicates that the Cu layer is completely intercalated. 

The effective intercalation of a thin Cu layer, and the formation of 
the graphene/Cu/Ir(lll) occurs at 550°C, is identified via strong 
modifications of the C Is and Ir 4/ emission lines [Fig. S2(a,b)]. 
Intercalation leads to a shift of the C Is peak to higher binding energy; 
it can be fitted with two components at 284.69 eV and 285.01 eV 
[Fig. S2(c)] [the integral intensity of the C Is line is restored to the 
value equivalent to the one for graphene/Ir( 1 1 1 ) ] . In case of the Ir 4/ 
spectra, the energy splitting between bulk (b) and interface (i) com- 
ponents is reduced from 537 meV to 463 meV, and the intensity of 
the interface component is strongly suppressed compared to that of 
the graphene/Ir(lll) system [Fig. S2(d)]. 

Although the STM data from the graphene with an intercalated Cu 
layer [Fig. l(b-d)] may appear similar to those from graphene/ 
Ir( 1 1 1 ) at first glance, there are clear differences in detail which reveal 
the effect of intercalation. The large scale STM image shown in (b) 
demonstrates the formation of a moire structure on top of Cu/ 
Ir(lll), which has the same periodicity of (10 X 10)graphene/(9 X 
9)Cu/Ir(lll). Hence we conclude that a pseudomorphic Cu layer on 
Ir( 1 1 1 ) is formed after intercalation. However, contrary to the results 
for graphene on Ir( 1 1 1 ) , a variation of the tunnelling voltage during 
STM imaging of graphene/Cu/Ir(lll) does not lead to an inversion 
of the imaging contrast - graphene is here always imaged in direct 
contrast [Fig. l(c,d)]. The two atomically-resolved STM images were 
acquired at (c) Ut = +0.3 V and (d) Ut = —0.3 V (upper part) and 
Ut — —0.9 V (lower part). In Fig. 1(d) the tunnelling voltage was 
changed during scanning "on-the-fly" showing the absence of any 
change of the contrast in the moire cell on an atomic scale. 

Figure 1(e) shows the structure of the graphene/Cu/Ir( 111) system 
as derived from DFT calculations. A (10 X 10) unit cell graphene 
lattice was placed on a (9 X 9) five layer Ir( 1 1 1 ) oriented slab, and the 
graphene and the topmost Ir layers were then allowed to relax. The 
resulting in-plane lattice constants of graphene are close to those of 
the free-standing species. For the structural evaluation of the inter- 
calated system, the topmost metal layer (S) in Fig. 1(e) was replaced 
by Cu. The Cu atoms were found to take up the positions of the Ir 
atoms, and the distance between the topmost metal layer and gra- 
phene in the ATOP position was found to be 3.581 A for graphene on 
Ir and 3 . 1 22 A for the intercalated Cu layer, a considerable reduction, 
signaling an enhanced interaction of carbon with the Cu atoms. 
Other structural parameters are given in Table Tl in the supplement- 
ary material. 

The conservation of long range periodicity of the system after 
intercalation of Cu is also shown by the LEED images of the gra- 
phene/Cu/Ir(lll) system [upper inset of Fig. 1(b)]. However, while 
the periodicity of the spots reflecting the existence of the graphene 
moire structure is similar to that for graphene/Ir(lll), the symmetry 
of spot intensities is reduced from sixfold to threefold. Hence we 
directly infer a lowering of the local symmetry, which means that 
the sublattice symmetry for the two carbon atoms in the graphene 
unit cell is broken. Thus we expect a strong modification of the 
electronic structure of graphene in the vicinity of E P . 

This modification is revealed in detail using ARPES and a 
comparison of the results with our DFT calculations. Such data 
demonstrate deviations from the band structure expected for 
free-standing or electronically decoupled graphene, and because of 
the sharp and well-separated Cu 3d band spectral features, permit to 
identify details of band hybridization. Consider the photoemission 
intensity maps of graphene/Ir(lll) and graphene/Cu/Ir(lll) in the 
vicinity of E P in Fig. 2. For graphene/Ir(lll) [Fig. 2(a)] the clear 
linear dispersion of the graphene n states is observed around E F with 
energy gaps at higher binding energies due to avoided crossings 
among the main n band and the replica bands which appear due 
to the additional periodicity of the moire lattice 21 . By extrapolation 
we extract a position of the Dirac point at about 100 ± 20 meV above 
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Figure 1 | STM images of graphene on Ir(lll) and pseudomorphic Cu/Ir(lll). (a) 3D view of the STM image of graphene/Ir(lll) (scanning 
parameters: 120 X 56 nm 2 , l/ r = +0.3 V, I T = 1.6 nA). Upper and lower insets show the LEED image (E p = 93 eV) and atomically resolved STM image 
(scanning parameters: 4.5 X 4.5 nm 2 , Ut = +0.5 V, I T = 10 nA) of this system, (b) 3D view of the STM image of graphene/Cu/Ir(lll) (scanning 
parameters: 175 X 82 nm 2 , Ut = +0.3 V, I T = 1.6 nA). Upper and lower insets show the LEED image (E p = 90eV) and an STM image of two graphene 
domains (scanning parameters: 60 X 24 nm 2 , U r = +0.3 V, I T = 1.6 nA). (c) and (d) present atomically resolved images of the graphene/Cu/Ir(lll) 
system obtained at Ut = +0.3 V and Ut = — 0.3 V(top)/ — 0.9 (bottom), respectively (scanning parameters: 10.1 X 10.1 nm 2 , I T = 1.6 nA for both 
images), (e) Top and side views of the graphene/metal(lll) structure. In case of Ir(lll) all layers consist of the same atoms, whereas for the Cu/Ir(lll) 
pseudomorphic system the metal (S) is the copper layer. In the top view, the circle, triangle, rhombus, and rectangular denote the ATOP, HCP, FCC, and 
BRIDGE positions for the carbon atoms, respectively. The big white rhombus marks the unit cell of the graphene/metal( 111) system. 



E F , corresponding to a slight p-doping of the graphene layer, in 
agreement with earlier ARPES data 21,31 . 

Intercalation of the Cu layer underneath graphene leads to a sig- 
nificant modification of the valence band states of the graphene layer 
[Fig. 2(b)]. Graphene on Cu/Ir(lll) now is n-doped with a position 
of the Dirac point at 0.688 eV below E P . Although the 7i band still has 
a linear dispersion near E p , a clear hybridization between the gra- 
phene n and Cu 3d valence band states in the 2-4 eV binding energy 
region is obvious. This manifests itself as a series of avoided crossing 
gaps between graphene- and Cu-d-derived valence band states [see 
also Figure 3]. Moreover, in the spectral function of the graphene n 
band, an energy gap of 0.36 eV appears at the Dirac point. This 
energy gap is more clearly resolved in photoemission data sets 
obtained at a photon energy hv = 40.81 eV presented in 
Fig. 2(c,d), where ARPES intensity maps around the K point are 
shown in (c) for the r — K (left panel) and perpendicular to T — 
K (right panel) directions as well as (d) corresponding energy cuts at 
the respective binding energies marked in the figure. It is surprising 
that the present values for the energy of the Dirac point E D and the 
gap are quite different from those for graphene on bulk Cu( 1 1 1 ) (E D 
— E p = —0.3 eV, gap = 0.25 eV 32 ), and for a single intercalated Cu 
layer in between graphene and Ni( 1 1 1 ) (E D — E F = — 0.3 eV, gap = 



0.18 eV 24 ). The different width of the gap for the different graphene/ 
Cu interface can be probably connected with the various periodicities 
of the corresponding moire structures as was shown in Ref. 33. 
Assuming the pseudomorphic growth of Cu layer at the interface 
in all cases the energy gap has to increase in the row graphene/Cu/ 
Ni(lll) graphene/Cu(lll) graphene/Cu/Ir(lll). 

Discussion 

The appearance of an energy gap around the K point in the spectrum 
of the graphene n states is one of the intriguing problems in graphene 
interface studies, and its existence or absence, for example in gra- 
phene on SiC(0001), has been extensively discussed 34 ' 35 . For single 
layer graphene on those metals that do not destroy the Dirac cone, i.e. 
the noble metals, Pt or Ir, gap openings have been observed in photo- 
emission, and this has been attributed to sublattice symmetry break- 
ing induced by a superstructure that couples the states at the K and K' 
points. As we will show through DFT calculations below, the dom- 
inating factor that determines the magnitude of the energy gap is the 
strength of hybridization of the graphene-derived n states with the 
valence band states of the substrate, which leads to the broken sym- 
metry of the two carbon sublattices as shown already by the LEED 
data in Fig. 1(b). 
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Figure 2 | Electronic structure of graphene/Ir(lll) and graphene/Cu/Ir(lll) measured by ARPES. (a,b) ARPES intensity maps for the graphene layer 
onlr(lll) and Cu/Ir( 111), respectively, acquired along the T — K direction of the BZ of graphene with a photon energy hv = 65 eV. The inset in (a) shows 
a graphene-derived BZ with the corresponding directions, (c) ARPES intensity maps for graphene/Cu/Ir( 111) obtained in the vicinity of the K point along 
r — K (left) and perpendicular to it (right) directions of BZ of graphene with photon energy hv = 40.81 eV. (d) The corresponding constant energy cuts 
extracted from the complete 3D ARPES data set for the graphene/Cu/Ir(lll) system. 



The analysis of structural features is shown in Fig. 1(e). Here the 
interface metal layer (S) is either Ir or Cu on Ir(lll). A pseudo- 
morphic arrangement of the Cu layer on Ir(lll) is assumed in the 
structure optimization procedure, a reasonable assumption in view 
of the LEED and STM data. At the energy minimum, a corrugation of 
the graphene layer on Cu/Ir(lll) of 0.229 A is found, which is much 
lower than the one of 0.307 A for graphene/Ir(l 1 1) 30 . The calculated 
carbon site projected partial density of states (PDOS) for the gra- 
phene-derived 7i states in graphene/Cu/Ir(lll) shows that graphene 
is n-doped with the Dirac point at 0.45 eV binding energy and a band 
gap of 0.15 eV [Fig. 4(a)], whereas the experiments give values of 
0.688 eV and 0.36 eV, respectively [Fig. 2(c)]. The computed DOS 



was obtained by broadening of the originally calculated data for the 
density of states obtained for the graphene/Cu/Ir(lll) supercell in 
the slab geometry: this broadening gives rise to a reduction of the 
energy gap, which would be 0.26 eV from the unbroadened data (see 
Fig. S4 from supplementary material). A comparison of the calcu- 
lated PDOSs for the Cu layer and the graphene layer indicates the 
existence of the hybridization between Cu 3d and graphene n states 
in the binding energy range of 1.5-3.5 eV [Fig. 4(a) and Fig. S5 from 
supplementary material] . 

The effect of hybridization can be visualized in real space in the 
charge distribution across the moire unit cell; the correspond- 
ing difference in electron density, Ap(r), defined as a difference of 
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Figure 3 | Analysis of the electronic structure of graphene/Cu/Ir(lll). (Upper panel) Comparison of the the experimental and calculated electronic 
structure of graphene/Cu/Ir(lll) and graphene/Cu(lll), respectively, in the energy and k- vector ranges where hybridization of graphene n and Cu 3d 
states is observed. (Lower panel) Orbital character decomposition of the valence band states of graphene in the vicinity of the Dirac point. See Fig. 1 (e) and 
Fig. S 1 in the supplementary material with the corresponding text for the definitions of the respective position of the carbon atoms. 



densities for graphene/Cu/Ir(lll) and those for the separate layers in 
the system, is shown in Fig. 4(b). The charge distribution picture for 
the FCC and HCP high symmetry sites is different from the one 
characteristic for graphene/Ir(lll) and similar to the situation for 
graphene/Ni(lll) or graphene/Rh(lll) where a strong interaction 
for these positions is observed. However, here the effect is weaker 
because the hybridization occurs between the valence states at higher 
binding energies. This charge distribution was then used to model the 
STM images. A simulated image corresponding to a bias voltage of Ut 



= —0.3 V is shown as an inset of Fig. 4, which is in good agreement 
with results presented in Fig. 1(d). As mentioned above, the ATOP 
positions of the graphene/Cu/Ir(lll) system are imaged as bright 
areas corresponding to direct imaging contrast; this situation prevails 
over a range from +0.3 V and —0.9 V [see Fig. S6 in the supple- 
mentary material]. This can be connected with the similar PDOS for 
different high-symmetry positions in the graphene moire structure for 
this system [Fig. 4(a)] and that in this case the real topography con- 
trast is prevailed in the STM images of graphene/ Cu/Ir(l 11). 
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Figure 4 | Results of the electronic structure DFT modelling for graphene/Cu/Ir( 111), (a) Carbon-atom projected site-resolved partial density of states 
for the graphene-derived n states. Greyed plot shows the corresponding DOS for the free-standing graphene. The band gap region is marked by the 
shadow area. Inset shows the simulated STM image corresponding to the experimental results obtained at [7 r = — 0.3 V. The high-symmetry positions are 
marked by the symbols similar to Fig. 1. (b) Side view of the graphene/Cu/Ir(lll) system with the corresponding difference electron density, Ap(r) = 
Pgr/Cu/Hm)(r) - P; r (iii)M _ Pc,,(ni)(r) - p gr (r), plotted in units of e/A 3 . 



Because of the narrow line widths of the Cu 3d-induced bands in 
photoemission, the graphene/Cu/Ir(lll) system is an excellent one 
to study band hybridization. Five features altogether can be distin- 
guished in the binding energy region from 1.8 eV to 4.0 eV. The 
avoided crossings of these bands with those derived from the gra- 
phene n band are clearly observed [shown enlarged in Fig. 3, upper 
panel, left], and the change of band character at around 2.15 eV, 
starting from mostly d-derived at k = —0.6 A -1 (with respect to 
the K point) to mostly graphene p z -derived is evident (topmost grey 
feature). Moreover, there are clear differences in the amount of 
hybridization, reflected in the magnitude of the avoided crossings. 
Because of the computational cost of deriving the band structure of 
the complete system (large unit cell of the moire structure), we 
restrict ourselves to an assignment of bands calculated for a Cu slab 
layer with the same lattice constant than graphene, i.e. 2.464 A. These 
data are shown next to the experimental bands in Fig. 3. Because of 
the slab model, a multitude of bands appear in the calculations; the 
main weight of the bands are indicated by color bands. Their relation 
to the different atomic d states from which they arise is indicated next 
to the calculation. The calculated bands are close to the experimental 
bands. 

Having shown the occurrence of hybridization of the graphene n 
and Cu 3d valence band states through photoemission, we dem- 
onstrate that this process is responsible for the appearance of the 
energy gap in the electronic structure of graphene at the Dirac point. 
Hybridization leads to the intermixing of the n orbitals with 3d 
orbitals of the corresponding orbital character (symmetry) depending 



on the carbon atom in the unit cell. According to the picture of 
Fig. 4(b), which shows the charge distribution difference, the most 
intensive interaction is observed around the FCC and HCP places. 
Here, one of the carbon atoms is placed above the interface Cu atom 
and hybridization between C" 9 p z and Cu 3d z i orbitals is observed. 
The second carbon atom is placed either above the hep or fee hollow 
site of Cu/Ir(lll), and here the hybridization between QS ccllcp p z and 
Cu 3d XZiyz orbitals occurs. For free-standing graphene the electronic 
states originating from two carbon sublattices are degenerate around 
the K point. The interaction of electronic states of carbon atoms with 
the 3d electronic states of different symmetry of the Cu layer leads to 
the lifting of this degeneracy and to the opening of the band gap at the 
K point. Symmetry breaking is also reflected in the photoemission 
data in the region of the Cu d bands. Along the T — K line, the group 
of the A:- vector in free-standing graphene is C 3v , such that carbon p z - 
derived bands and those from Cu d xzyz would not be allowed to 
interact. However, an avoided crossing between these states, indi- 
cative of symmetry reduction, is apparent in the data of Fig. 2. The 
symmetry must therefore be further reduced, on account of the lattice 
mismatch leading to the moire structure. 

This conclusion is supported by analyzing the weight of the con- 
tributions from different orbitals to the bands below and above the 
band gap as a function of energy and electron wave vector. The lower 

panel of Fig. 3 shows the weight of the p z '° p + d^"^ (top) and the 
pj" + d x "( y f (bottom) hybridized states, represented by the mag- 
nitude of the symbols. It is clear that the p z + d z i hybrid state 
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dominate the bands below the band gap, while the p z + d XZtyz one 
dominates the bands above. This clearly shows how hybridization 
leads to sublattice breaking evidenced, for example, by the reduced 
symmetry of the LEED pattern in Figure 1 . Comparing the graphene/ 
Cu/Ir(lll) system with the original graphene/Ir(lll) interfaces, the 
interaction or hybridization between graphene valence band state 
and Ir d orbitals is weaker compared to the Cu d orbitals and hence 
the sublattice symmetry breaking is more pronounced in the latter 
case. This to our knowledge is the first case where such a detailed 
assignment of orbitals contributing to the states around the band gap 
has been achieved. Similar analysis for graphene on other substrates 
in which gap opening occurs will permit an understanding of the 
interplay of substrate and graphene states, and will shed light on the 
as yet elusive correlation between bonding strength, magnitude of 
band gap, and its consequences for the structural arrangement. 

In summary, STM and photoemission data from an intercalated 
layer in between graphene and Ir(lll), interpreted on the basis of 
state-of-the-art DFT calculations, permit a detailed analysis of the 
changes in doping and band gap opening in the graphene n bands 
upon intercalation, in terms of state-specific hybridization between 
the different Cu d and graphene n bands. The strong spatial and 
energy overlap of the valence band states leads to the lifting of the 
sublattice symmetry of the carbon atoms which manifests itself as an 
opening of the band gap at the Dirac point. Analyses of this kind 
permit to predict the arrangement of the electronic states of graphene 
in other graphene/metal interfaces. 

Methods 

Preparation of graphene/Ir(lll) and graphene/Cu/Ir(lll). The graphene/Ir(lll) 
system was prepared in ultrahigh vacuum system for ARPES experiments according 
to the recipe described in details in Refs. 27, 30, 36 via cracking of propylene: T = 
1 100 C, p — 5 X 10~ 8 mbar, t — 30 min. This procedure leads to the single-domain 
graphene layer on Ir( 1 1 1 ) of very high quality that was verified by means of LEED and 
STM (in the separate experiments). Intercalation of the 1 ML-thick Cu layer and 
formation of the graphene/Cu/Ir{lll) intercalation-like system was achieved via 
annealing of the graphene/Ir( 111) sample with the thin pre-deposited copper layer on 
top. The process of intercalation was monitored by measuring of the "live" C Is and Ir 
4f spectra of the system and formation of the graphene/Cu/Ir(l 1 1) system was 
detected at 550 C. The quality and cleanness of the system in ARPES experiment was 
verified by LEED and XPS/ ARPES, respectively. The base vacuum was better than 8 X 
10~ n mbar during all experiments. All ARPES and STM measurements were 
performed at room temperature. 

ARPES experiments. The ARPES measurements at photon energies of hv — 65 eV 
and hv — 94 eV were performed at the BESSY II storage ring (Berlin) in the 
photoemission station using PHOIBOS 100 2D-CCD hemispherical analyzer from 
SPECS. Experiments with He II radiation (hv = 40.81 eV) were performed in the 
SPECS demo lab using a FlexPS system with PHOIBOS 150 2D-CCD analyser. In 
both cases a 5-axis motorized manipulator was used, allowing for a precise alignment 
of the sample in k space. The sample was azimuthally pre-aligned in such way that the 
tilt (BESSY II) or polar (FlexPS) scans were performed along the T — K direction of 
the graphene- derived BZ with the photoemission intensity on the channelplate 
images acquired along the direction perpendicular to F — K. The final 3D data set of 
the photoemission intensity as a function of kinetic energy and two emission angles, 
KEfon, angle^ angle 2 ), were then carefully analyzed. 

STM experiments. The STM measurements were performed in constant current 
(CC) mode. In this case the topography of the sample, z{x,y), is studied with the 
corresponding signal, tunneling current (I T ), used as an input for the feedback loop. 
The STM images were collected with SPM Aarhus 150 equipped with the 
KolibriSensor™ from SPECS 30,37 with a Nanonis Control system. In all 
measurements, a sharp W-tip was used which was cleaned in situ via Ar + -sputtering. 
In the STM images the tunnelling bias voltage, U T , is referenced to the sample and the 
tunnelling current, I Ti is collected by the tip, which is virtually grounded. 

DFT calculations. The crystallographic model of graphene/Ir(lll) and graphene/ 
Cu/Ir(l 1 1) presented in Fig. 1 was used in the DFT calculations, which were carried 
out using the projector augmented plane wave method 38 , a plane wave basis set with a 
maximum kinetic energy of 400 eV and the PBE exchange- correlation potential 39 , as 
implemented in the VASP program 40 . The long-range van der Waals interactions 
were accounted for by means of DFT-D2 approach proposed by Grimme 41 . The 
studied system is modelled using supercell, which has a (9 X 9) lateral periodicity and 
contains one layer of ( 10 X 10) graphene on a five-layer slab of metal atoms. Metallic 
slab replicas are separated by ca. 20 A in the surface normal direction. To avoid 



interactions between periodic images of the slab, a dipole correction is applied 42 . The 
surface Brillouin zone is sampled with a 7 X 7 X 1 fc-point mesh centered the F point. 
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